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ABSTRACT 

Starting as highly relativistic collimated jets, gamma-ray burst outflows gradually decellerate and become 
non-relativistic spherical blast waves. Although detailed analytical solutions describing the afterglow emission 
received by an on-axis observer during both the early and late phases of the outflow evolution exist, a calculation 
of the received flux during the intermediate phase and for an off-axis observer requires either a more simplified 
analytical model or direct numerical simulations of the outflow dynamics. In this paper we present light curves 
for off-axis observers covering the long-term evolution of the blast wave, calculated from a high resolution 
two-dimensional relativistic hydrodynamics simulation using a synchrotron radiation model. We compare our 
results to earlier analytical work and calculate the consequence of the observer angle with respect to the jet 
axis both for the detection of orphan afterglows and for jet break fits to the observational data. We find that 
observable jet breaks can be delayed for up to several weeks for off-axis observers, potentially leading to 
overestimation of the beaming corrected total energy. When using our off-axis light curves to create synthetic 
Swift X-ray data, we find that jet breaks are likely to remain hidden in the data. We also confirm earlier results 
in the literature finding that only a very small number of local Type Ibc supernovae can harbor an orphan 
afterglow. 

Subject headings: gamma-rays: bursts - hydrodynamics - methods: numerical - relativity 



1. INTRODUCTION 

According to the standard fireball shock model, gamma- 
ray burst (GRB) afterglows are the result of the interac- 
tion between a decelerating relativistic jet and the surround- 
ing medium. Synchrotron radiation is produced by shock- 
accelerated electrons interacting with a shock-generated mag- 
netic field. The radiation will peak at progressively longer 
wavelengths and the observed light curve will change shape 
whenever the observed frequency crosses into different spec- 
tral regimes, when the flow becomes non-relativistic and 
when lateral spreading of the initially strongly collimated out- 
flow becomes significant (see e.g. Zhang & Meszaros 200^ 
[Piran 2005; Meszaros 2006 for recent reviews). 

Analytical models have greatly enhanced our understanding 
of GRB afterglows. Such models rely on a number of sim- 
plifications of the fluid properties and radiation mechanisms 
involved. Both at early relativistic and late non-relativistic 
stages spherical symmetry can be assumed. At first lat- 
eral spreading of the jet has not yet set in and the beam- 
ing is so strong that a collimated outflow is still observa- 
tionally indistinguishable from a spherical flow. Eventually 
the outflow really has become approximately spherical. Self- 
similar solutions for a strong explosion can be applied, the 
Blandford-McKee (BM, Blandford & McKee 1976) solution 
in th e relativistic regime and t he Sedov-Taylor-von Neumann 
(ST. I Sedovlll 9591: iTavloril 19501) solution in the non-relativistic 
regime. However, in order to include lateral spreading of 
the jet and to calculate the light curve for an observer not 
located on the axis of the jet, the downstream fluid profile 
has usually been ap proximated by a homogeneous slab (e.g. 
Rhoads.,19 99: Kumar & Panaitesculi2000: Granot et aL..2002: 
Waxman| ]2004: Oren et al.l l2004l) . Structured jet models ex- 



ist where the effect of lateral expansion is estimated and fluid 
quantities like Lorentz factor and density de pend on the angle 
of the flow with respect to the jet axis (e.g. IRossi et al.ll2008l 



also lGran"o3l2007 l and references therein). 

However, to gain an understanding of afterglow light 
curves during all stages of jet evolution and for off-axis 
observers, large scale multi-dimensional simulations are re- 
quired. Over the past ten years various groups have combined 
one-dimensional relativistic hydrodynamics (RHD) simula- 
tions with a radiation calculation (e. g. Kobavashi et al. 199a 
'Dow nes et all l2002t iMimica et all 120091: I Van Eerten et all 
2010ai) . Thanks to specialized techniques such as adaptive- 
mesh-refinement (AMR), jet simulations i n more than one di- 
mension have also be come fea sible (e.g. iGranot et al.l l 2001l : 
IZhang & MacFadveri 2006; M eliani et alj|2007h 

In this work we present the results of a high-resolution two- 
dimensional RHD simulation covering the full transition from 
relativistic to non-relativistic flow. We use the simulation re- 
sults from iZhang & MacFadvenll2009l (hereafter ZM09), but 
we now calculate for the first time detailed light curves for 
off-axis observ ers. Th e simulations have been performed with 
the RAM code (IZhang & MacFadyen 2006). 

Off-axis observations of GRB afterglows are observation- 
ally relevant for a variety of reasons. Following the first obser- 
vations of GRB afterglows, it was immediately realized that, 
if the afterglow emission was less strongly beamed than the 
prompt emission, many orphan afterglows should in principle 
still be observable even when the prompt emission was not 
visible because the observer^as positioned too far away from 
the jet axis ( Rhoads 1997. 1 19991) . Detailed light curves from 
simulations help constrain the expected rate of occurrence of 
orphan afterglows and help determine to what extent orphan 
afterglows can possibly remain hidden in observations of type 
Ibc supernovae. This in turn helps to constrain the fraction of 
type Ibc supernovae that can in principle be linked to GRBs. 

For observers that are not too far off-axis but are still lo- 
cated within the jet opening angle, the prompt emission still 
remains visible and they will observe an afterglow light curve 
that shows a jet break when the jet edges become visible, an 
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effect which is enhanced when lateral spreading becomes sig- 
nificant. The shape of this jet break is expected to depend 
strongly on the observer angle. If the observer is not on- 
axis, each edge of the jet becomes visible at different times 
and the corresponding break is split in two, or at least be- 
comes smoother. This effect may acco unt for the difficulty 
in detecting a jet break for many GRBs jRacusin et alj|2009l : 
lEvans et al.ll2009l) . 

This paper is structured as follows. In section lZTI we briefly 
review the RHD simulation setup and methods we have ap- 
plied in Z M09 and that form the basis for this paper as well. 
In section |Z21 we describe how the radiation is calculated for 
an observer at an arbitrary angle with respect to the jet axis. 
The resulting light curves are presented in section|3] They are 
put into context first by a comparison to light curves calcu- 
lated from the relativistic BM solution for on-axis observers, 
followed by a comparison to light curves from a simplified ho- 
mogeneous slab model for off-axis observers. We then apply 
the simulation to two different observational issues. In sec- 
tion|4]we use our computed light curves to generate synthetic 
Swift data, to which we then fit broken and single power laws 
in order to probe the extent to which X-ray jet breaks can be 
hidden in the SWIFT data due to off -axis observer angle. In 
section |5] we confirm the result from iSoderberg 61111(120061) 
that only a very small number of local Ibc supernovae can pos- 
sibly harbor an orphan GRB afterglow, now using our simu- 
lation instead of a simplified analytical model for comparison 
with the observations. We discuss the results presented in this 
paper in section |6] The mathematical details of the analytical 
model with which we have compared our simulation results 
are summarized in the Appendix. 

2. METHODS 

2.1. Hydrodynamic Model 

We have used the two-dimensional RHD simulation first 
presented in ZM09 as the basis for our calculations. This 
simulation was p erformed using the RAM a daptive mesh re- 
finement code (Zhang & MacFadyen 2006). Ram employs 
the fifth-order weighted essentially non-oscillatory (WENO) 
scheme (Jiana & Shu 1996) and uses the PARA MESH AMR 
tools (.MacNeice et al...2000.) from FLASH 2.3 (.FrvxeU et all 
l2000h . 

The simulation takes a conic section of the Blandford- 
McKee (BM) analytical solution (Blandford & McKee 1976) 
as the initial condition, starting from a fluid Lorentz factor 
directly behind the shock front equal to 20. The isotropic en- 
ergy of the explosion was set at Eiso = 10^^ erg. The jet half 
opening angle 9q = 0.2 rad (11 .5°), leading to a total energy in 
the twin jets of Ej w 2.0 x lO''' erg. The circumburst proton 
number density is taken to be homogeneous and set at « = 1 
cm"-*. The pressure p of the surrounding medium is set at a 
very low value compared to the density p (p= 10~^^pc^, with 
c the speed of light) and will therefore not be dynamically 
important. Under these conditions, the starting radius of the 
blast wave is equal to Rq « 3.8 x lO'^ cm. 

A spherical grid (r,d) was used with < r < 1.1 x 10'^ 
cm and Q < 6 < tt/I. At first, 16 levels of refinement are 
used, with the finest cell having a size of Ar w 5.6 x 10'^ cm 
and A9 ss 9.6 x 10"^ rad. The maximum refinement level is 
gradually decreased (but always kept at least at 1 1) during the 
simulation, making use of the fact that the blast wave widens 
proportionally to in lab frame time to keep the number of 



cells radially resolving the blast wave approximately constant. 

The most important dynamical results from ZM09 are the 
following. They find that very little sideways expansion takes 
place for the ultrarelativistic material near the forward shock, 
while the mildly relativistic and Newtonian jet material fur- 
ther downstream undergoes more sideways expansion. When 
taking a fixed fraction of the total energy contained within an 
opening angle as a measure of the jet collimation it is found 
that sideways expansion is logarithmic (and not expo nential, 
as used by some early analytic models such as that of Rhoadsl 
1999). This sideways expansion sets in approximately at a 
time tg calculated from plugging 7 = l/9j into the BM so- 
lution, where 7 the fluid Lorentz factor directly behind the 
forward shock and 9j the original jet opening angle. For the 
simulation settings described above, tg ss 373 days, measured 
in the frame of the burster The jet becomes nonrelativistic 
and the BM solution breaks down at f ^ f^g w 970 days. The 
time Ink is estimated by equating the isotropic equivalent en- 
ergy in the jet to the rest mass energy of the material swept up 
by a spherical explosion, assuming the jet moves at approxi- 
mately the speed of light. The transition to spherical flow was 
found to be a slow process and was found from the simulation 
to take until St^R to complete. After that time the outflow can 
be described by the Newtonian Sedov-von Neumann-Taylor 
(ST) solution. 

2.2. Calculation of Off-Axis Afterglow Emission 

A large number of data dumps (2800) from the hydrody- 
namic simulation are stored. These data dumps are then used 
to calculate the synchrotron emission from all the individual 
fluid elements at the lab frame time of each data dump. A 
single emission time corresponds to different arrival times at 
the observer for different parts of the fluid due to light travel 
time differences. The radiation from the data dumps is binned 
over a number of observer times. The size of each fluid el- 
ement is given by AV = r^sin6'ArA0A0 in 3D. In the 2D 
axisymmetric simulation the flow is independent of cf) and, if 
the observer is positioned on-axis, the (p symmetry allows for 
considering just the 2D elements AV = Inr^ sin 6ArA9 pro- 
vided by the data dumps. When calculating emission for off- 
axis observers, however, fully 3D data must be created from 
the 2D fluid elements by extending the data in the direction. 
The fluid elements are split into smaller elements along the 
direction of angular size A^ to account for differences in 
relativistic beaming and Doppler shifts of emission observed 
from different angles. In practice we have started with angular 
resolution such that rsin0A0 is comparable to Ar and rA9, 
taking into account that the arrival time between the close and 
far edge of the fluid element in the direction should stay very 
small. We found that the light curves are not very sensitive to 
the resolution in the direction, even when it is decreased 
tenfold. 

The sy n chrotr on emission itself is calculated following 
ISari et al.l ([T998h . We sum over the contributions of the in- 
dividual fluid elements. In the frame comoving with the fluid 
element, the spectral power peaks at 

, /16\^ p— 1 ajmpC^ , , 

P' = 0.88 — —^n'B' (1) 

V 3 y 3p-l %nq, 

where p is the slope of the power law accelerated electron dis- 
tribution, aj is the Thomson cross section, is the electron 
mass, c is the speed of light, is the electron charge, «' is 
the comoving number density and B' is the comoving mag- 
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netic field strength. The field strength B' is determined from 
the comoving internal energy density e\ using B' = y^ege^Svr. 
Here is a free parameter that determines how much energy 
is converted into magnetic energy near and behind the shock 
front. The shape of the spectrum is determined by the syn- 
chrotron critical frequency v'^ and the cooling frequency v',.. 
These frequencies are set according to 



3 

16 



qeB' 

p—l n'nieC^ J nieC ' 



(2) 



where e^. parameterizes the fraction of the internal energy den- 
sity in the shock-accelerated electrons, and 
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The cooling break is therefore estimated by using the duration 
of the explosion as a measure of the cooling time: t is the lab- 
frame time and 7 is the Lorentz factor of the fluid element in 
the lab frame (i.e. the frame of the unshocked medium) used 
to translate to the time comoving with the fluid element. The 
emitted power at a given frequency depends on the position of 
that frequency in the spectrum with respect to f,', and v'^ and 
the relative position of i',', and v'^ (fast cooling when < j/,', 
and slow cooling when < u'^). The spectrum consists of 
connected power law regimes. For example, for a fluid el- 
ement the emitted power at observer frequency i/ iv' in the 
comoving frame) between synchrotron break and cooling in 
the slow cooling case is given by P'(v' lv'„,)^^~'''*l'^ . The c om- 
plete shape of the spectrum is given in lSari et"an (1199 8*) and 
can also be found in the Appendix. The power in the observer 
frame is then obtained by applying the appropriate beaming 
factors and Doppler shifts to power and frequency. Finally 
the received flux is calculated by taking into account the lumi- 
nosity distance and redshift (the latter is also used to translate 
between observer frequency and comoving frequency). 

In our calculations we have set p = 2.5, e^, = = 0.1. These 
are typical values found for afterglows. Different redshifts 
have been calculated, but in this paper we only present re- 
sults where we have ignored redshift (i.e. z = 0) and set the 
observer luminosity distance at di = 10^^ cm. If the redshift 
is increased the features of the light curves stretch out to later 
observer times. Light curves computed for a range of observer 
redshifts will be presented in an upcoming publication. 

The main radiation results from ZM09 for an on-axis ob- 
server are the following. The jet break due to lateral expan- 
sion was found to be weaker than analytically argued, while 
the jet break due to jet edges becoming visible was stronger 
than expected from the simplest analytical models, although 
not unexpected from calculations taking limb-brightening into 
account. The weaker jet break due to lateral spreading can 
be understood from the lateral spreading being logarithmic 
instead of exponential as has been often assumed in analyt- 
ical models (see ZM09, Fig. 3). The long transition time 
to the nonrelativistic regime for the blast wave was already 
mentioned in the previous section. When the blast wave has 
become nonrelativistic the counterjet is no longer beamed 
away from the observer It becomes distinctly visible around 
tcj = 2(1 +z)tNR, with the ratio between flux from counter and 
forward jet peaking at 6 at 1 GHz at 3800 days for the simu- 
lation settings (at z = 1). 

3. NUMERICAL RESULTS 



We first present the main results for both small and large 
observer angles, looking both at the light curves and the cor- 
responding temporal slopes. An overview of light curves is 
shown in Fig. [T] where we have plotted multi-frequency light 
curves spanning from 10"^ to 10'^ Hz. The temporal slopes 
for the lowest frequency 10^ Hz and the highest frequency 
10^^ Hz are separately plotted in Figs. |2]and[3] For the on- 
axis results we estimate the jet break, a combination of lateral 
spreading and jet edges becoming visible, to occur around 3.5 
days (the 7 days mentioned in ZM09 is for z = 1). Direct 
comparisons at the same frequency of different observer an- 
gles are shown in figures |6] and [TT] for 8.46 GHz. As the 
observer angle increases, the jet break splits into two for ob- 
servers still inside the jet. Once the observer is positioned at 
the jet edge only one break remains that is significantly post- 
poned. This effect is similar across all frequencies. The steep- 
est drop in slope is the one associated with the edge of the jet 
furthest from the observer and can therefore be estimated to 
occur around 



f; = 3.5(1 +z)£, 
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where £,.50.53 is the isotropic equivalent energy in units of 
10^^ erg, «! is density of the medium in units of cm"-', 6'o is 
the jet half opening angle, and d„bs is the observer angle rel- 
ative to the jet axis. Jet breaks can be used to estimate the 
opening angles of GRB jets. It is usually assumed in GRB 
afterglow modeling that the observer is on the jet axis. How- 
ever, if the observer is near the edge of the jet, the jet opening 
can be overestimated by a factor of up to 2, and the beaming- 
corrected total energy can be overestimated by a factor of up 
to 4. For a typical observer at dots ~ 2(?o/3, the beaming- 
corrected energy can be overestimated by a factor of ^ 3. The 
observational implications of this effect will be further dis- 
cussed in Section |6] 

At high observer angles, the rise of the light curve is post- 
poned until the point where relativistic beaming has weakened 
sufficiently for the observer to be in the light cone of the ra- 
diating fluid. Due to limb-brightening the drop in temporal 
slope following a jet break initially overshoots its asymptotic 
value. After that it starts to change again due to the onset 
of the transition into the nonrelativistic regime and the rise 
of flux from the counterjet, before it finally settles into its 
asymptotic value for the nonrelativistic regime. In order to 
put the simulation results in context and to differentiate be- 
tween the break due to lateral spreading and the break due to 
the edges becoming visible we will compare the simulation 
results against the BM solution for a hard edged jet without 
lateral spreading in the next subsection. 

3.1. Afterglow emission for on-axis observer - 
hydrodynamic simulation versus analytic model 

In Fig. Hlwe show a comparison between on-axis spectra 
at 1 day in observer time, calculated from the simulation and 
from an analytical description u sing the BM solution plus syn- 
chrotron emission (IVan Eerten & Wiiers 2009). The observed 
time is well before the jet break and before significant lateral 
spreading or slowing down of the jet has occurred. With the 
dynamics for both the simulation and the BM solution still 
being nearly equal, the figure therefore mainly shows the dif- 
ference between the two approaches to synchrotron radiation. 
The differences below the cooling break Vc are marginal and 
can be attributed to the absolute scaling of the emitted power. 
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Fig. 1 . — Simulated light curves for various observer angles and frequencies. On the top row we have small observer angles, 0.0, 0.1 and 0.2 radians from left 
to right. On the bottom row we have large observer angles, 0.4, 0.8 and 1.57 radians from left to right. Observer frequencies from lO' up to lO'^ are plotted. 
For each frequency and angle a curve is plotted both with and without the contribution from the counterjet. The green vertical lines in the top plots indicate jet 
break time estimates using equation |4] The legend in the bottom right plot refers to all plots. Large observer angles are truncated at earlier time to show only 
observation times completely covered by the simulation. 
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Fig. 2. — Temporal decay index a assuming F oc f The observer fre- 
quency is lO' Hz. The vertical lines at 3.5 days indicate the jet break time 
estimates for observers at 0, 0.1 and 0.2 radians (from left to right), using 
equation |4] 



Fig. 3. — Temporal decay index a assuming f oc f The observer fre- 
quency is lO" Hz. The vertical lines at 3.5 days indicate the jet break time 
estimates for observers at 0, 0. 1 and 0.2 radians (from left to right), using 
equation |4] 



The difference beyond is significantly larger. The reason 
for this is that th e simulation foll ows the approach to elec- 
tron cooling from lSari et alj (Il998h . where the cooling time is 
globally estimated by setting it equal to the duration of the 
explosion, whereas .Van Eerten & Wiiers (2009.) build upon 



'Granot & Sari' ('2002") and calculate the local cooling time for 
each fluid element, which is given by the time passed since 
the fluid element has crossed the front of the shock. When 
the cooling time is calculated locally, the transition between 
pre- and post-cooUng is smooth, with areas of the fluid further 
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Fig. 4. — direct comparison between on-axis spectrum from simulation 
(solid curve) and BM solution with locally calculated cooling times (dashed 
curve), omitting self-absorption. The spectra are taken at 1 day in observer 
time, so well before the jet break. The leftmost vertical dotted line denotes 
the analytically calculated position of v„,, the rightmost that of Vc (calculated 
for the BM solution). 



tion of this break and the jet break due to lateral spreading. 
The exact light curves shown in Fig. |5]do not cover the entire 
observer time span because the BM solution ceases to be valid 
around a shock Lorentz factor of 7 ^ 2 (and slightly earlier for 
a local cooling calculation). The light curves are truncated at 
the point where this would start to affect the observed emis- 
sion. 

From the figure it can be seen that including lateral spread- 
ing has the effect that the jet break becomes steeper and starts 
slightly earlier (which confirms the low resolution compari- 
son shown in Van Eerten et al. 2010b). The figure also shows 
the strong overshoot in steepening of the light curve following 
the jet break , also seen in Fig. [3] This overshoot has also been 
discussed in lGranot etall (2001) and ZM09. 

The difference between the detailed local treatment of elec- 
tron cooling and the global treatment of cooling is important 
when applying simulation results to actual data: it should 
be kept in mind that the simulation light curves systemati- 
cally underestimate the flux beyond the cooling break. Elec- 
tron cooling aside, the long term qualitative behavior of the 
light curve is fully captured by the simulation, with the re- 
sults covering not only the relativistic regime but also the non- 
relativistic regime and the transition in between. 



downstream making the transition before those at the front. 
Having the cooling time set globally results in a sudden global 
transition between the pre- and post-cooling regimes instead, 
with the same asymptotic spectral slope but a different value 
for the cooling break frequency and therefore a different abso- 
lute scaling of the flux for v > Vc- Note that even for a global 
cooling time, the sudden transition in the emission frame will 
still get smeared out in the observer frame. 
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Fig . 5 . — Direct comparison of simulation results and heuristic description 
based on the Blandford-McKee exact solution. Simulated light curves are 
shown for observer frequencies lO'^ and lO'^ Hz (top and bottom solid line 
respectively). Analytical light curves are drawn for global (dashed) and local 
{dotted) coohng respectively. The top line of each pair is for a spherical model 
whereas the lower line is for a conic section only. The 3.5 days estimate for 
the jet break is indicated with a vertical dotted hne. 

Figure |5] shows a direct comparison between the on-axis 
light curves obtained from the simulation and light curves 
from the same analytical model as before. If no lateral spread- 
ing is assumed, a conic section of the spherically symmetric 
BM solution can be used to show the difference between the 
jet break due to the edges becoming visible and the combina- 



3.2. Afterglow emission for off-axis observer - 
hydrodynamic simulation versus analytic model 

No exact solution exists that fully includes lateral spreading 
of the jet. In the Appendix we describe a simplified analytical 
model that approximates the behavior of the jet and allows 
us to calculate the observed flux for an observer at an arbi- 
trary angle . Many such models exist in the literature ( see e.g . 
Oren et al. (2004); WaxmanI (12004 : ISoderberg et"an (|2006l): 
Huang et al. (2007) etc.) and our model is not strongly dif- 
ferent. Its distinguishing features are that it smoothly con- 
nects the relativistic BM solution to the nonrelativistic ST so- 
lution and that a conservative approach to lateral spreading is 
used where the jets start to spread at the speed of sound (and 
therefore logarithmically) upon approaching the nonrelativis- 
tic regime. 
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Fig. 6. — Direct comparison between simulation results (solid lines) and 
analytical model (dashed lines) for different opening angles, for radio fre- 
quency 8.46 Ghz. 

In figure |6] we show a comparison between off-axis light 
curves generated using the analytical model and light curves 
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Fig. 7. — Light curves at 8.46 Ghz for an observer at Oots = 90°. The sim- 
ulation and simplified model curves are repeated from fig. |6] If we ignore the 
fluid velocity in the lateral direction for the puipose of the radiation calcula- 
tion, we get the slightly lower curve labeled no lateral beaming. The emission 
is still well in excess of a light curve from a radiation calculation that takes the 
BM profile instead of the simulated fluid profile as input, labeled Blandford- 
McKee. Only by completely omitting the emission contribution from matter 
that has spread sideways out of the original jet opening angle, we get a flux 
level that is initially comparable to that of the exact BM solution and the 
simplified model. This is shown by the truncated cone curve. 

calculated from the simulation. Qualitatively both simulation 
and model show the same features. Quantitatively however, 
the differences are substantial. The simulation light curves 
peak earlier than the model light curves and do so at lower 
peak luminosity. At early times, the emission from off-axis 
simulation light curves is higher than that from the corre- 
sponding model light curves. The early time slopes for sim- 
ulation and model light curves are similar However, the fur- 
ther off-axis the observer, the later the observer time at which 
the simulation provides full coverage. In the figure, we have 
truncated the light curves at the observer time before which 
radiation from the blast wave with 7 > 20 would have been re- 
quired. For two blast wave jets viewed sideways {dots = 90°), 
for example, full coverage starts beyond 100 days. This 
means that, even though the initial slopes agree between simu- 
lation and model, the early time shape of off-axis light curves 
in reality will be largely dictated by the initial shape of the 
blast wave, which does not need to be anything like the BM 
solution. Collapsar jet simulations indicate the existence of a 
cocoo n around the emerging je t (see e .g., Zhang et al. 2003, 
lIOOliMorsonvetal. 2007; Mi zuta & Al ov 2009). 

We can understand why the off-axis light curves from the 
simulation are initially brighter than those from the simpli- 
fied model by looking at one of the angles in more details. In 
Fig. |7]we have again plotted the simulation and model light 
curves for an observer at 0„bs = 90° (1.57 rad), together with 
a number of variations. We have now also included a light 
curve where we continue using the BM solution to determine 
the local fluid conditions, instead of the dynamical simulation 
results, but otherwise proceed as if we were reading the fluid 
quantities from disc (because of this, the curve also serves as 
a consistency check on the radiation calculation itself). The 
same approach has been used to generate the BM global cool- 
ing light curve in Fig. |5] The curve initially lies significantly 
below the simulation light curve. It also lies above the sim- 
plified model curve. The flux level of this BM curve is deter- 
mined in part by the numerical resolution that we assume. In 



the plot we have used a resolution similar to that for the sim- 
ulation curve, which initially resolves the radial profile with 
approximately 17 cells. Increasing the resolution moves the 
BM curve closer to the simplified slab curve, and not to the 
simulation curve. The difference between the BM curve and 
the simulation curve is real and we have added to the plot two 
hybrid simulation / model curves to make clear the cause of 
this difference. 

First, when we completely ignore the velocity vg in the an- 
gular direction for the purpose of calculating the emission but 
otherwise still use the simulation dynamics, we find that the 
resulting light curve, labeled no lateral beaming in the figure, 
initially lies somewhat below the simulation curve before the 
two eventually merge. This tells us that part of the observed 
flux level is caused by beaming towards the observer of mate- 
rial spreading sideways, but that this is not the main cause of 
the difference between simulation and the hard edge jet mod- 
els. At late times beaming no longer plays any role and the 
two curves are indeed no longer expected to be different. 

The main reason for the difference is shown by the second 
additional curve. When calculating the light curve labeled 
truncated cone in the figure we have omitted the contribu- 
tion to the radiation of any material that has spread sideways 
outside the original jet opening angle. The resulting curve 
lies very close to the BM light curve at first, before becom- 
ing orders of magnitude lower than all other curves. The late 
time behavior is as expected, for then only a small fraction 
of the energy and particle density is still contained within the 
original opening angle; the actual simulation flow has become 
roughly spherical. The early time behavior and the similarity 
between the truncated cone and BM light curves is more rel- 
evant. It demonstrates that the light curve for an off-axis ob- 
server is dominated by the emission from material that has 
spread sideways out of the original jet opening angle even 
though the energy of the material is very little and the side- 
ways spreading is not yet dynamically important. In hind- 
sight, the fact that the material on the side of the jet dominates 
the observed radio flux can easily be understood. It is not so 
much due to the fact that it moves a little faster towards the ob- 
server, for the Vg component to the beaming is not that strong 
(as we have shown above). It is instead due to the fact that 
the radial velocity component vv drops quickly outside of the 
original jet opening angle and as a result the material outside 
the original jet opening angle is not beamed away from the 
observer as much as the material in the original jet cone. By 
contrast, for an on-axis observer the opposite is true and for 
a long time the received flux is dominated by emission from 
material inside the origina l jet opening angle. Thi s has been 
demonstrated explicitly by I Van Eerten et al] (1201 Oh) . 

In the above sections we did not discuss the effects of syn- 
chrotron self-absorption. We will postpone addressing syn- 
chrotron self-absorption, which is not included in the simula- 
tion, until section |5j 

4. APPLICATION: HIDDEN JET BREAKS? 

A large number of X-ray afterglow light curves have been 
obtained by the Swift satellite since it was launched in 2004 
(IGehrels etalJl20 07i. In a surprisingly large number of cases, 
these light cu rves fail to sh ow a clearly discernable jet break 
dRacusin et al. 2009; Eva ns et al.ll2009h . Using our simula- 
tion results as a basis to generate synthetic Swift data sets for 
observers positioned at different angles from the jet axis, we 
show that the effect of observer position on the temporal evo- 
lution inferred from the data can be profound and sometimes 
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render the jet break difficult to detect. 

4.1. Procedure for creating synthetic data 

The synthetic data sets that we produce should be com- 
parable to those produced by the on-line Swift repository 
dEvans et al.l2007h . Also they should have data points at a suf- 
ficiently late time that, if this were actual data, the jet break 
would be considered missing and not merely delayed. We 
therefore make sure that we have data up to at least 10 days, 
in accordance with the c riteria for their 'complete' sample set 
by iRacusin et al.l (120091) . The observed on-axis jet break for 
our simulation occurs roughly around three days. Synthetic 
hg ht curves hav e also been created from an underlying model 
bv lCurran et al.l (2008) (who find that even broken power law 
models observed on-axis can occasionally be mistaken for a 
single power law decline) and we follow the same procedure 
as described in that paper, changing only the time span and 
adding an additional late time data point if necessary. We 
then have: 

• Constant counts and 1 a fractional error of 0.25 per data 
point. 

• 94 minute orbits (47 min on/off due to Swift's low-Earth 
orbit). 

• Fractional exposure drops from 1.0 to 0.1 after one day 
(when Swift is usually no longer dedicated completely 
to observing the burst). 

• Rate cut off at 5 x 10""* cts/s. 

• Observed number of cts/s is scaled to 0. 1 at 1 .0 day. 



10" 



simulation 
synthetic Swift data 



simuiation 
syntlietic Swift data 




observer time (s) 



Fig. 8. — Synthetic Swift data generated from simulation liglit curve for an 
on-axis observer and p = 2.5. Tlie simulation flux has been scaled to a count 
rate of 0. 1 cts/s at 1 day observer time. 

We start generating data points from 3 x 10"* s and continue 
until the rate drops below 5 x lO"'* cts/s. This starting point 
is chosen such that we have full coverage from t he simulation 
at all o bserver angles under consideration. Like iCurran et al.l 
(l2008h . we increase the number of counts per bi n to a num- 
ber we ll in excess of the numbers mentioned by IE vans et aTl 
( 120071) . This has no physical significance but is used to gen- 
erate synthetic light curves containing around 30 data points 
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10" 



10" 



observer time (s) 



Fig. 9. — Synthetic Swift data generated from simulation light curve for an 
on-edge observer (i.e. at 0.2 rad) and p = 2.5. The simulation flux has been 
scaled to a count rate of 0. 1 cts/s at 1 day observer time. 



(The synthetic curves from Curran et al. contain more data 
points because they use an earlier starting time). If needed a 
late time data point is added to ensure that at least one data 
point is observed after ten days. The last one or two data 
points, where the count rates are less than lO""* cts/s get a 
larger fractional error of 0.5. Out of the thousands of syn- 
thetic curves that have been generated, two randomly selected 
example synthetic light curves are shown in figures [8] and |9] 

We generate light curves for observer angles 0.00, 0.02, 
0.04, 0.06, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18 and 0.2 radians. 
Note that by fixing the count rate at 0. 1 cts/s after one day, we 
end up comparing on-axis observations to off-axis observa- 
tions that are relatively brighter (i.e. corresponding to closer 
GRBs). As can be seen in Fig. [T] off-axis light curves are less 
bright than on-axis curves for the same physical parameters. 

4.2. Fitting procedure and results 

We follow ICurran et al.1 ( 120081) again when fitting the 
synthetic data automatically using the simulated annealing 
method to minimise the of the residuals. The data are first 
fit to a single power law, then to a sharply broken power law. 



F(t)=N 



(t/tb)-'" if t<ti„ 
(t/tb)-°''- iff >?/,. 



(5) 



Curran et al. use a smooth po wer law, but a sharply broken 
power law is also used by Racu sin et al.l (120091) . After each fit 
the count rates of the data points are re-perturbed from their 
original on-model values and a Monte Carlo analysis using 
1000 trials is used to obtain average values and la Gaussian 
deviations of the best fit parameters and F-test probabilities. 
This process is repeated for the list of observer angles men- 
tioned previously. We have set p = 2.5. 

The F-test is a measure of the probability Fprob that the de- 
crease in associated with the addition of the two extra pa- 
rameters of the broken power law, and f/,, is by chance 
or not. When Fprob ^ 10"^ a single power law is commonly 
favored, when 10"^ < Fpi-„b < 10"^ neither is favored but a sin- 
gle power law is usually presumed as the simpler model and 
when Fpiob ^ 10"^ a broken power law is favored. 

The results of the fitting procedure are shown in tables [Tl|2] 
and fig. [TO] In table [T] any slope a and its error Aa means 
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observer angle (rad) 

Fig. 10. — Results for F-test as function of observer angle. Note that a 
typical observer angle is 0.13. 



dobs 


a 




CK2 


tb (10^ s) 


0.00 


1.69 ±0.040 


0.81 ±0.16 


3.32 ±0.79 


2.3 ±0.80 


0.02 


1.68 ±0.037 


0.76±0.17 


3.14±0.73 


2.1 ±0.83 


0.04 


1.63 ±0.038 


0.81±0.15 


2.96 ±0.61 


2.1 ±0.72 


0.06 


1.78 ±0.036 


0.84 ±0.16 


3.05 ±0.58 


2.1 ±0.70 


0.08 


1.52 ±0.040 


0.78±0.19 


2.26 ±0.37 


1.4 ±0.54 


0.10 


1.51 ±0.040 


0.80 ±0.20 


2.07 ±0.26 


1.2 ±0.45 


0.12 


1.49 ±0.042 


0.91 ±0.26 


1.95 ±0.29 


1.2 ±0.55 


0.14 


1.65 ±0.035 


1.07±0.16 


2.36 ±0.52 


2.0±1.1 


0.16 


1.46 ±0.042 


1.00 ±0.28 


1.95 ±0.47 


1.6 ±0.89 


0.18 


1.40 ±0.040 


0.91 ±0.36 


1.85±0.50 


1.6 ±0.97 


0.20 


1.42 ±0.036 


0.86±0.21 


2.01 ±0.45 


1.9±1.0 



TABLE 1 

Average temporal power law slopes and jet break times fj, 
FOR 1000 Monte Carlo iterations per observer angle So^j. 
Here a denotes the single power law slope, ai the pre-break 
broken power law slope and q2 the post-break broken power 
law slope. 



dobs 


log(F,,„i,) 


single 


ambiguous 


broken 


0.00 


-7.1 ±1.7 





130 


870 


0.02 


-6.9 ±1.6 





148 


852 


0.04 


-6.6 ±1.6 





202 


798 


0.06 


-7.2 ±1.6 





60 


940 


0.08 


-5.2 ±1.4 


6 


510 


484 


0.10 


-4.5 ± 1.4 


22 


698 


280 


0.12 


-3.6 ± 1.4 


129 


767 


104 


0.14 


-4.6± 1.5 


21 


641 


338 


0.16 


-2.8 ± 1.3 


274 


675 


51 


0.18 


-3.1 ±1.4 


259 


686 


55 


0.20 


^.4 ±1.5 


37 


663 


300 



TABLE 2 

Average log(F,„.ofc) values and classifications based on the 
criteria described in the text. for each angle, the 
classifications add up to 1000 Monte Carlo runs. 



that if an observer fits data to a swift dataset with similar 
count rate and duration and if the swift data was produced 
by an explosion that is accurately described by our numerical 
model, then that observer would find a slope that lies within 
Aa of a (the errors therefore do not reflect on the accuracy 
of the Monte Carlo run, which has converged to far greater 
precision). The analytically expected slopes for an on-axis 



observer and p = 2.5 are 1.375 before and 2.125 after the jet 
break. These values are not even reproduced within 1 a, which 
for the post-break slope can largely be attributed to the over- 
shoot in slope directly after the break (see also section[3]i. Al- 
though the inaccuracy of the pre-break slope inferred from 
(synthetic) light curves will be smaller when earlier time data 
is available, which is often the case for S wift light curves, a 
system atic difference will remain (see also'Johan nesson et al.l 
120061) . This emphasizes the importance of numerical model- 
ing for the proper interpretation of Swift data. The situation 
gets worse when we move the observer off-axis. As Fig. [TO] 
shows, even though a jet break is clearly detected on-axis for 
our physics settings, it can become hard to distinguish from 
a single power law for an observer positioned at Oj/2. The 
average observer angle one would expect when observing jets 
oriented randomly in the sky is 29j/3 (for small jet opening 
angles, and assuming that the observer angle lies between zero 
and the jet opening angle). Larger observer angles will lead to 
orphan afterglows that we discuss separately below. It there- 
fore follows that a significant number of jet breaks may re- 
main hidden in the data due to the jets not being observed di- 
rectly on-axis. In table |2]the classifications for the individual 
Monte Carlo iterations are also counted. Assuming again that 
jets are oriented randomly in the sky and that jets are only ob- 
served up to observer angles equal to dj, we can calculate how 
often an afterglow with the physical parameters of the simu- 
lation would be classified as showing a jet break. For each 
angle 6*, we have classified out of 1000 synthetic curves 
as showing a jet break. This means that we will classify the 
afterglow described by the simulation as showing a jet break 
only X;,n/ sin(6l,)/ X),- 1000sin(6',) x 100% « 29% of the time. 
This value is only a very rough estimate, for it depends on the 
Fprob ^ 10"^ criterion, that is to some extent arbitrary. Also, 
as noted before, the off-axis curves are relatively brighter due 
to the fixed count rate. The requirement of having data up to 
ten days introduces a selection effect as well. 

For all results above, it should be kept in mind that they 
have been obtained for a single half opening angle of 0.2 ra- 
dians (approx. 11.5°) which is relatively large (although not 
extremely so and within the range of jet opening angles ob- 
servationally inferred from Swift data). This results in a later 
jet break than a smaller opening angle would lead to. On 
the other hand we have set z = 0, which again moves the jet 
break to earlier observer times and thereby compensates for 
the large opening angle. The general effect of higher observer 
angles is that both jet edges become observable at different 
times. The reason that a broken power law did not always 
produce a significantly better fit than a single power law is 
mainly because the full drop in count rate associated with the 
further edge got pushed out beyond ten days for off-axis ob- 
servers (Section O. For strongly collimated jets with small 
opening angles this effect is therefore expected to be less se- 
vere. 

The often large difference between the physical parameters 
like p (affecting the slope of the light curve) and 6j (affect- 
ing the break time) used in a model and the values for these 
parameters when inferred from synthetic data created from 
that model has been discussed in detail bv lJohannesson et"an 
(,2006.) . They also include the observer angle as a model pa- 
rameter, but do not discuss it further in their paper. 

5. APPLICATION: ORPHAN AFTERGLOW SEARCHES 

The existence of orphan afterglows is an important and 
general prediction of current afterglow theories. Regard- 
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less of the GRB launching mechanism and the initial baryon 
content of the jet, eventually synchrotron emission from a 
decelerating baryonic blast wave should be observable for 
any observer angle. For this reason, various groups have 
looked for orphan afterglows, both at the optical and radio 
frequencies (e.g . Levinson et al, 2002' Gal- Y am et al. 2006; 
ISoderberg etal.l 2006; Malacr ino et allliOOTl) . Few positive 
detections have been reported and surveys and archival stud- 
ies have mainly served to establish con s traints on GRB rates 
and beaming factors. Soderberg e t al.l (120061) . for example, 
conclude from late time radio observations of 68 local type 
Ibc supernovae (SNe) that less than ^ 10% of such SNe are 
associated with GRBs, and constrain the GRB beaming fac- 
tor to be ^(1 -cos^y)"') < 10"^. A lower limit to the beaming 
factor of ^(1 -cos6'y)~') > 13 is provided by iLevinson et al.l 
(l2002i) . 




10^ 10^ 10^ lO'* 

observer time (days) 



Fig. 1 1 . — VLA late time radio limits (3cr) for 66 local type Ibc of super- 
novae compared against simulation results. All supernova redshifts have been 
ignored (the largest redshift, that of SN 1991D, is ^ 0.04). The fluxes have 
been rescaled to luminosities. All VLA observations were done at 8.46 GHz, 
the afterglow light curve is calculated at the same frequency. As in the rest of 
this paper, the simulation jet half opening angle is 11.5°. 

Such estimates require a model describing the shape of off- 
axis light-curves, and are therefore sensitive to model assump- 
tions. Eventually, comparing observations and detailed simu- 
lations like the one described in this paper will place the most 
accurate observational limits on orphan afterglow characteris- 
tics. A large number of simulations is required to fully explore 
the afterglow parameter space. We can however, use the sin- 
gle simulation of this paper that has typical v alues for the ex- 
plosion parameters to confirm the result from lSoderberg et al] 
12006) that their sample of 68 supernovae observations are aU 
significantly fainter than a standard afterglow viewed off-axis. 
This confirmation is shown in fig. (TT] where we have plot- 
ted 66 supernovae radio upper limits (omitting SN 1984L and 
SN1954A, which were not observed at 8.46 GHz, from the 
original 68) together with our off-axis simulated light curves. 
Note that the jet half opening angle in our simulation is 1 1 .5°, 
whereas Soderberg et al. (2006) use 5°. The fact that the early 
time flux received by an off-axis observer is actua lly s tronger 
than analytically expected (as shown in section 13.21 where 
model and simulation are compared directly) only strength- 
ens the case made by Soderberg et al. 

A possible caveat to the above is that our simulation 




10^ 10^ 10^ lo" 

observer time (days) 



Fig. 12. — Analytically calculated light curves for different observer an- 
gles, with and without self-absorption. The observer frequency is set at 8.46 
GHz. Self-absorption influences the hght curves only until a few hundred 
days at most. At 8.46 GHz, the light curve for an observer at 90 degrees is 
the same with and without self-absoiption enabled. The jet half opening angle 
is 11.5°. The VLA late time radio limits for the local Ibc SNe are included 
as well. 

light curves do not include the effect of synchrotron self- 
absorption. Although we cannot completely rule out that this 
plays a role without actually calculating it, we can neverthe- 
less look at the effect of self-absorption on the model light 
curves, having already established that model and simulation 
lead to at least qualitatively similar light curves in section [J!2l 
In Fig. [12] we show model light curves with and without syn- 
chrotron self-absorption, calculated as explained in the ap- 
pendix. The figure shows that the effect of self-absorption is 
initially significant for an on-axis observer but becomes less 
pronounced for observers further off-axis. For an observer at 
90°, the light curves with and without self-absorption are ef- 
fectively identical. Aside from the minimal differences due 
to the analytical model assumptions, the main difference be- 
tween this figure and fig. 1 from Soderberg et al. is due to 
the different jet opening angles. For our wider jet opening an- 
gle, only the two earliest supernovae lie clearly above the 90° 
curve. 

6. SUMMARY AND DISCUSSION 

In this paper we present broadband GRB afterglow light 
curves calculated assuming synchrotron emission from a 
high-resolution relativistic jet simul ation in 2D. We have ex- 
panded the work presented in Zhan g & MacFadyeiil (120091) to 
include observers positioned off the jet symmetry axis, both 
at small and large angles. For the jet simulation we have 
used the RAM adaptive-mesh-refinement code, starting from 
the Blandford-McKee analytical solution and letting the jet 
evolve until it has reached the Sedov-Taylor stage and has 
decoUimated into a nearly spherical outflow. We have im- 
pleme nted synchrotron radiation as described in iSari et all 

jTHi. 

When put in the context of analytical light curve estimates, 
our simulations show the following: 

• For an on-axis observer, the jet break from a 2D simula- 
tion including lateral spreading of the jet is seen earlier 
than that of a hard-edged jet. However, the break due to 
the jet edges becoming visible still dominates the shape 
of the light curve. 
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• We compared a description of electron cooling that 
takes into account the local cooling time since the 
shocked electrons passed the shock front to a descrip- 
tion that uses a global cooling time estimate (as we have 
done in the simulation). The latter approach underesti- 
mates the observed cooling break frequency and there- 
fore the post-break flux as well. 

• The simulation light curves show that simplified homo- 
geneous slab analytical models are qualitatively correct 
as long as they include a clear contribution from the 
counterjet for observers at all angles. 

• Contrary to what has thus far been assumed in simpli- 
fied analytical models and even though lateral spread- 
ing of the jet is initially not dynamically important, the 
received flux for an off-axis observer is strongly dom- 
inated by emission from material that has spread later- 
ally outside the original jet opening angle. This is due 
to the fact that material outside the original jet cone has 
slowed done considerably in the radial direction and is 
therefore not beamed away from the observer as much 
as material closer to the jet axis. 

• Moving the observer off-axis splits the jet break in two. 
The steep drop in slope only occurs after the break asso- 
ciated with the farthest edge and this break can thus be 
postponed until several weeks after the burst, even for 
observers positioned still within the jet cone. The late 
break time can be estimated by using the sum of the 
observer angle and the jet half opening angle instead of 
just the jet half opening angle. 

In addition to these direct numerical results, we have pre- 
sented two applications of our numerical work, one with ob- 
servers positioned at small and moderate observer angles and 
one with observers at large observer angles. 

• X-ray Light curves for observers at small observer an- 
gles are re levant for satellites such as Swift. Rece nt au- 
thors (e.g. iRacusin et al ."2009*: "E vans et al.ll2009l) have 
noted a lack of jet breaks visible in the data. In order 
to check whether the observer angle can cause the jet 
break to remain hidden in the data we have performed a 
Monte Carlo analysis where we created synthetic Swift 
data out of simulation light curves for different observer 
angles. Observational biases and a Gaussian observa- 
tional error were included. Broken and single power 
laws were then fit to the synthetic data sets. We found 
that it is not difficult to bury a jet break in the data for an 
off-axis observer For our explosion parameters, even 
an observer at an angle of Oj/2 will not find a signifi- 
cantly better fit for a broken power law than for a single 
power law, leading to a missing jet break. For a random 
observer angle within the cone of the jet, a synthetic 
light curve created from our simulation will only show 
a discernable jet break 29 % of the time. Although our 
simulation is somewhat atypical in its above average 
initial opening angle of the jet, these results neverthe- 
less imply that the observer angle has a strong influence 
on the interpretation of X-ray data. This holds even 
for observers still within the cone of the jet, observer 
angles that have usually been ignored and considered 
practically on-axis. 



• As a second app l ication we have confirmed the result of 
ISoderberg et al.l (120061) that a sample of 68 nearby type 
Ibc SNe cannot harbor an off-axis GRB, at least for the 
typical afterglow parameters that we have used for our 
simulation. This confirms the observational restrictions 
placed on orphan afterglow rate and beaming factor by 
these and other authors. The fact that at early times the 
light curves for off-axis observers are actually brighter 
than anal ytically expected, even strengthens the conclu- 
sions of lSoderberg et al.1 (120061) . 

Recent deep late-time optical observations by iDai et al] 
(|2008) have detected jet breaks in several bursts and they 
suggest that the lack of jet breaks in Swift bursts is due to 
the lack of well sampled light curves at late times. There- 
fore they conclude that the collimated outflow model GRBs 
is still valid. However, the non-detection of a jet bre ak or a 
break a t very late times in several GR Bs (050904, Cenk o et al. 
'2010b; 070125, Cenko et al."2010b; 080319B, Cenko etal. 
2010b; 080721, Starling et al. 2009; 080916C, Greineret al 
2009; 090902B. iPandev et al.i 12010 : 090926A. iCenko et al. 
,2010a) seems to infer a huge amou nt of released energ y 
(> 10^2 erg) in fliese bursts. This leads ICenko et all (l2010bllah 
to propose a class of hyper-ener^etic GRBs and challenge the 
magnetar model for GRBs (Usov 1992; Duncan & Thompso^ 
1992; Thompson et al. 2004; Uzdenskv & MacFadven :20()7F 
.Komissarov & Barkov ^007; Bucci antini et al. 2009) . Al- 
though hyper-energetic GRBs might exist, we propose an al- 
ternate explanation, in which the observer is simply off-axis. 
Note that an observer is more likely to be off-axis than on- 
axis. A typical observer sees the burst from Oohs ~ "^Sd/li so is 
closer to the jet edge than jet axis. Applying an on-axis model 
to a jet that is seen off-axis will overestimate the total energy 
release by a factor of up to 4 because the jet break can be de- 
layed by a factor of up to ^ 6 (see Eq. |4|i. We emphasize that 
when the jet opening angle is corrected for off-axis observer 
angle, the inferred energy release of those hyper-energetic 
events can be revised downwards by factors of several lessen- 
ing tension with magnetar models for the GRB central engine. 

These are the main conclusions of the work presented in 
this paper. In addition to this, the current results also raise 
a number of issues that need to be addressed in future work. 
Our simulation results can be generalized by performing ad- 
ditional simulations. And given the quantitative differences 
between simulation and simplified analytical models (such as 
the homogeneous slab model described in the Appendix), this 
is expected to result in different constraints on orphan after- 
glow rate and beaming factor We have so far ignored self- 
absorption when calculating light curves. A simplified homo- 
geneous slab approximation indicates that self-absorption is 
not expected to play a large role for observers at high angles. 
The applicability of our analytical model is however limited, 
especially in view of our finding that the early flux received 
by an off-axis observer is dominated by emission from mate- 
rial that has spread sideways and has slowed down more than 
material on the jet axis. It is conceivable that emission from 
this material has different spectral properties. Finally, the sig- 
nificant difference in observed flux between two approaches 
to synchrotron cooling that we have discussed emphasize the 
importance of a detailed model for the microphysics and radi- 
ation mechanisms involved. 

We thank Peter A. Curran for allowing us the use of 
his computer code for synthesizing and fitting Swift data 
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APPENDIX 

SUMMARY OF ANALYTICAL MODEL 

Here we provide a short summary of our analytical model. In this model, light curves are calculated by numerical integration 
over an infinitesimally thin homogeneous blast wave front. The received flux is given by 

^ I'd'r ,,f-' (Al) 



ATidlJ 72(1-/3^)2 

when we ignore the redshift z- Here ej^, is the comoving frame emissivity and /i the angle between observer direction and local 
velocity. The dependence of the beaming and emissivity on the observer time tots is kept implicit (see also eqn. IA8t the volume 
integral needs to be taken across different emission times). Assuming the radiation is produced by an infinitesimally thin shell 
with width IS.R we get 

For every observer time we integrate over jet angles 9 and 4> (we define 6 such that it is not the angle between observer and local 
fluid velocity, but that between local fluid velocity and jet axis), while taking into account that radiation from different emission 
angles. 

The emissivity can be calculated from the local fluid conditions, which we know in turn in terms of emission time /«. For the 
blast wave radius we have by definition: 

R = cjp,hdte, (A3) 

Where the subscript sh indicates shock velocity. From the shock jump conditions it follows for arbitrary strong shocks that 

e',„ = (j-l)n'mpc\ (A4) 

The comoving downstream number density n' in both the relativistic and nonrelativistic regime is given by 

n' = 4no7, (A5) 

with 7 — >^ 1 in the nonrelativistic limit. We assume this equation to remain valid in the intermediate regime as well. This is not 
implied by the expression above, where we have kept implicit the dependence on the fluid adiabatic index (which changes from 
4/3 to 5/3 over the course of the blast wave evolution). 

We set the width of the shell at a single emission time by demanding that the shell contains all swept-up particles, leading to: 

4^ [Rf(tf)-Rh(tf)] R\ = 4/3nR\o [Rf(tf)-Rh(tf)] = (A6) 

iZ7 

where we have used n = jn' = AriQj^, again assumed valid throughout the entire evolution of the fluid. The subscript / denotes 
the front of the shock and the subscript b denotes the back of the shock. Setting the shock width through the number oiF particles 
is to some extent an arbitrary choice, and we could also have used the total energy which would have yielded a different width 
(since the downstream energy density profile is different from the downstream number density profile). The width of the shell 
AR in equation lA2l has to take into account the emission time difference between the front and back of the shell and is given by 



^R=\Rf{tf)-Rb{tb)\ = \Rf(tf)-Rf{tt) 



1- ' 



127(ffc)2 



(A7) 



Because the shell is very thin, Rjiti,) w Rf(tf)- /SshcAt. We integrate over emission arriving at a single observer time, and for 
given values of /i and tabs we have 

tohs = tf-^lRfitf)/c = tt-^lRb{tb)/c, (A8) 
which yields IS.R = t^tcj [i when differentiated. Combining the above, we eventually find 



A7?= i (A9) 



For the shock velocity we have 

Wshlsh)BM=i ^ ' "° I = CBMf;^/2. (/3,,,7,,)sr = --1.15- ^) ■tf/'=CsT-tf^', (AlO) 

in the BM and ST regime respectively. We artificially combine the two simply by adding them (after squaring): 

Phi = ClM'e' + Clrt;''' • (All) 
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Note that the BM quantities depend on Eis,,, while the ST quantities depend on Ej. The two are related via Ej = Eiso0j/2. Here Ej 
is the total energy in both jets, and 6j the half opening angle of a jet. The fluid Lorentz factor in the relativistic regime is related 
to the shock Lorentz factor via 7^ = 7j/,/2, while the fluid velocity in the non-relativistic regime is related to the shock velocity 
via f3 = 3/4/3 5/,. We therefore construct a relationship between emission time and fluid velocity similar to that between emission 
time and shock velocity: 

= \cl^t;^ + ^cjjtfi' . (A12) 

We assume that the jet does not spread sideways throughout the relativistic phase of its evolution. However, at some point the 
blast wave must become spherical. For if we keep the opening angle fixed, but do take Ej to dictate the ST solution instead of 
Eis„ we would underestimate the final flux by integrating over an integration domain that is too small. We will assume that the jet 
starts spreading sideways when it has reached the nonrelativistic phase and we take this moment to be given by 

{h)BM=^^tt,R = 2'I^Cy^. (A13) 
At this point, the jet starts spreading sideways with the speed of sound c,, leading to 

de 

R-7— = Cs. (A14) 

In the nonrelativistic regime dtots and dt^ are identical. The ST solution for the speed of sound for adiabatic index 5 /3 is given 
by Cj = r/ y/l^t, leading to 

61 = 61 , + -^ In — , (A15) 
■' %/20 tNR 



until spherical symmetry is reached. By contrast, Rhoads '99 take Vl ss TT{6j + Cst' /ctef' (where f' time in the comoving frame, 
and Vl a solid angle) as the starting point. 
From the local fluid conditions the local emissivity can be calculated. In the case of slow cooling, we define 

1/3 



(l-p)/2 



, X (l-p)/2 / , X -pI2 



v' < v'. (A16) 

The definition for fast cooling is analogous (see also lSari et al.lll998l) . The peak emissivity is given by 

, p-1 V3ql , , 

e,„ ^ ■ (^17) 

2 nieC 

The synchrotron break frequency i',', depends on its corresponding critical electron Lorentz factor 7^, leading to 

= Z (7m) ^ > 7m = 7 2' ^^^^^ 

For the cooling break frequency an identical relation between frequency and electron Lorentz factor holds. The critical Lorentz 
factor is now given by 

67rme7C 

which follows from the electron kinetic equation when synchrotron losses dominate over adiabatic expansion and the cooling 
time is approximated by the lab frame time since the explosion. 

Synchrotron self-absorption is included in the model using the assumption that emission and absorption occur in a homoge- 
neous shell. The solution to the linear equation of radiative transfer then dictates that we need to replace eq. IA2l bv 



7; =^7^, (A19) 



^ [ d0d(l)R^sine—il-e-^), (A20) 
rdl J Uy 



F = 

A-Kd] 

The optical depth r w -aiyA/jQ Emissivity and absorption translate between frames using e'^, = — /3fi)^ei, and a'^,/ = 
ai//7(l respectively. In our simplified model we calculate the absorption coefficient a'^, under the assumption that electron 
cooling does not influence it. This assumption is justified when the self-absorption break frequency I'a hes well below the cooling 

> This is an approximation that does not take into account that not all rays increasing the optical depth does not alter our finding that self-absorption 

cross the homogeneous slab along the radial direction. However, significantly ^^^^ pj^^ ^ ^^^^ YLA light curves generated from this model. 
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break frequency Vc, which is the case for all applications of the model in this paper. Also approximating the synchrotron spectral 
shape by just two sharply connected power laws we then find for the self-absorption coefficient: 

where k= 1/3 if < v'^ and K = -p/2 otherwise. 

Numerically speaking the integration procedure is as follows. First we tabulate R(te) for a given set of physical parameters, so 
that we do not need to estimate it analytically but can use its exact dependence on the fluid Lorentz factor instead. We integrate 
over 9 before we integrate over (j), and for each 9, (f> the angle between observer and fluid element is given by 

^ = sin 6* cos (/) sin 6*0^ + cos6'cos 9obs ■ (A22) 

Having tabulated R(te), we tabulate fi(te,R(te)) as well for a given value of tabs- Since is a monotonically increasing 
function of t^, we can unambiguously determine te(fi) from this table. When determining the value of the integrand at a given 
value of 9, </>, we can now calculate the local fluid conditions and emissivity via te{iJi{9, 0)). 
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